Spin-Polarized Transport in Ferromagnet-Marginal Fermi Liquid Systems 
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Spin-polarized transport through a marginal Fermi liquid (MFL) which is connected to two non- 
collinear ferromagnets via tunnel junctions is discussed in terms of the nonequilibrium Green function 
approach. It is found that the current-voltage characteristics deviate obviously from the ohmic be- 
havior, and the tunnel current increases slightly with temperature, in contrast to those of the system 
with a Fermi liquid. The tunnel magnetoresistance (TMR) is observed to decay exponentially with 
increasing the bias voltage, and to decrease slowly with increasing temperature. With increasing 
the coupling constant of the MFL, the current is shown to increase linearly, while the TMR is found 
to decay slowly. The spin-valve effect is observed. 

PACS numbers: 71.10.Ay, 75.25.+Z, 73.40.Gk 



Spin-polarized transport in magnetic hybrid nanos- 
tructures has been an active subject under investiga- 
tion in last decades, which is mainly motivated by po- 
tential applications in information technology. A new 
field coined as spintronics is thus emerging (for review, 
see e.g. Refs. Q, S S S 0)- The well-known char- 
acter in spintronics is that the current flowing through 
the structures depends sensitively on the relative orien- 
tation of the magnetization directions due to the spin- 
dependent scattering of conduction electrons. Among 
others, the magnetic tunnel junction (MTJ) is an im- 
portant family of spintronic devices^, 0- For these 
structures, Julliere [2j was the first to observe the tun- 
nel magnetoresistance (TMR) of 14% in Fe-Ge-Co junc- 
tions at 4.2 K. In 1995, Moodera et al. [9( made a 
breakthrough by observing reproducibly a large TMR as 
high as 24% at 4.2A and 11% at 295A'. Recently, clear 
spin-valve signals at 4.2A as well as at room tempera- 
ture have been observed in ferromagnet-normal metal- 
ferromagnet (FM-N-FM) all-metal structures |l(i||. Ear- 
lier theories on the spin-dependent transport in FM-N- 
FM j unctions E3 are based on the Fermi liquid theory, 
where interactions between electrons in the normal metal 
are treated on a mean-field level. There has been recent 
studies on the spin transport in FM-Luttinger liquid- 
FM tunnel junctions where the interactions between elec- 
trons are taken into account, and applied directly to 
carbon nanotubes |l2l but they are primarily aimed 
at one-dimensional interacting quantum wires. Besides, 
spin-polarized transport through an interacting quantum 
dot that is described by the Anderson model has also 
gained much attention[l4|. On the other hand, there ap- 
pear intriguing experimental and theoretical works on 
the spin-polarized transport in FM-high Tc supercon- 
ductor tunnel junctions recently (e.g. Refs.fl5|h It is 
thought that the anomalous normal state properties of 
high Tc cuprates in the optimally doped regime can be 
well described by the marginal Fermi liquid (MFL) 0, 
where the interactions between electrons in the cuprates 
are phenomenologically included in a one-particle self- 
energy due to exchange of charge and spin fluctuations. 



Therefore, the study on the spin-dependent transport in 
FM-MFL-FM tunnel junctions would be interesting, as it 
would be useful for understanding the transport proper- 
ties of FM-high Tc cuprate junctions in the normal state. 

In this paper, by using Keldysh's nonequilibrium 
Green function formalism, the spin-dependent transport 
in FM-MFL-FM tunnel junctions is investigated. It is 
observed that the current-voltage characteristics in this 
spintronic structure show non-ohmic behaviors, and the 
tunnel current increases slowly with temperature, which 
are in contrast to those of the structure with a Fermi liq- 
uid, showing that the interactions between electrons in 
the normal metal have remarkable effects on the trans- 
port properties. The TMR is found to decay exponen- 
tially with increasing the magnitude of bias voltage, and 
to decrease slowly with increasing temperature. With in- 
creasing the coupling constant A of the MFL, the current 
is shown to increase linearly, while the TMR is seen to 
decay slowly, implying that the interactions of electrons 
tend to suppress the TMR. In addition, the spin-valve 
effect is observed. 

Let us consider a MTJ in which the two FM electrodes, 
connected with the bias voltage V/2 and — V/2, respec- 
tively, are separated by a normal metal which is described 
by the MFL, as schematically depicted in Fig. 1. The 
molecular field in the left (L) FM is assumed to be 
parallel to the z axis, while the molecular field h/? in the 
right (R) FM is parallel to the z 1 axis which deviates the 
z axis by a relative angle 0. Tk aq (a = L, R) stand for the 
elements of the tunneling matrix between the a electrode 
and the central region. The tunnel current flows along 
the x axis and perpendicular to the junction plane. In 
the central region, the interactions between conduction 
electrons are supposed to be described phenomenologi- 
cally by a retarded one-particle self-energy due to the 
exchange of charge and spin fluctuations jig: 

£( E ) = A[ £ ln-|-*|z], (1) 

where x — max(|e|, UbT), E c is a cut-off energy and A is 
a coupling constant. When A = 0, the MFL junction re- 
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Green function in the central region. By applying Lan- 
grenth theorem |l7j | and Fourier transform, one may ob- 
tain formally 
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FIG. 1: (Color online) Schematic illustration of the double 
tunnel junction consisting of two ferromagnets (FM) and a 
marginal Fermi liquid (MFL) separated by insulating films, 
where T^ aq (a = L, R) stand for the elements of coupling 
matrix between the a electrode and the central region, and 
both magnetizations are aligned by a relative angle 6. 



covers the conventional Fermi liquid. For simplicity, the 
spin-orbital coupling in the MFL will be ignored. It is 
worthy of noting that the exact Hamiltonian of the MFL 
is not yet available. However, since the single-particle 
Green function is explicitly written down, the concrete 
form of the microscopic Hamiltonian is irrelevant. In the 
calculations, we just need to adopt a formal Hamiltonian 
such as a Fermi liquid with the electron operators under- 
stood as those of quasi-particles. Because the final results 
are all expressed by Green functions, we only need to use 
the MFL Green functions to replace the quasi-particle 
Green functions. 

By means of the nonequilibrium Green function, the 
tunnel current through the left electrode can be obtained 
by 



I L (V) = e(N L ) 



kqcr 



< 
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where Nl is the occupation number of electrons in the left 
electrode, G< a , Lkcr (t, t') = i(a] to .(i') £ V(*)) is the lesser 
Green function, a ka and Cko- are annihilation operators of 
electrons with momentum k and spin a (= ±1) in the left 
and central region, respectively. In order to get the lesser 
Green function, we define a time-ordered Green function 
G q ^ Lkff (M') = -i(T{at a (t') C(l<7 ,(t)}). In terms of the 
equation of motion, we have 



where G r , /(T (e) is the Fourier transform of the re- 
tarded Green function of electrons in the MFL of the 
central region, and G^ a , , a (e) is the corresponding lesser 

Green function, g^Lai 6 ) ana ^ 9kLa( £ ) are tne l esser 
and advanced Green functions for the uncoupled elec- 
trons in the left electrode. By defining r Q (e) q / - qcr ' = 
27rD(e)Tk Qq Tk Qq '5 (T(T ' with D(e) the density of states 
(DOS) in the a electrode and using the Fourier trans- 
form, after a tedious but direct derivation, Eq. J3} can 
be rewritten as 



x[f L (e)(G r (e)~G a (e)) + 



f cr Ml ) 
G < (e)}}, (4) 



where f a (s) is the Fermi function of the a electrode, and 
Tr is the trace over the momentum and spin space. Note 
that in Eq. 0} all Green functions, G r,a,< (e), are for elec- 
trons in the MFL of the central region, where G r ' a (e) are 
known with the presumed self-energy E(e) in the MFL 
[Eq. Q], say, G r (e) = [e - e k - ££ - E(e) + ir,]-\ 
where Eq, E(e) denote the coupling of MFL to the two 
ferromagnets and the retarded self-energy of the MFL, 
respectively, while G < (s) is unknown and needs to be 
obtained. 

To get the lesser Green function G < (e) of the central 
region, we invoke Ng's ansatz[l8j|: £< = E^£>, where 
££(e) = i[fL(e)T L (e + f- + <jM l ) + f R (e)RT R (s - 
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the rotation matrix, and M R = 
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5Mb^--r/2- Under this presumption, one may find even- 
tually that Eq. 0} becomes 
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where g^ La {t) = (ihj^ _ e kLcr ) with e kLa = e L (k) - 
(eV/2) — ctMl, £f,(k) the single-particle dispersion in the 
left electrode and Ml = g/J B IiL/2 (g- Lande factor, /i B : 
Bohr magneton), and G qo ./ q / CT (t — t') is the time-ordered 



The TMR ratio can be defined according to the current 
as usual El: 



TMR 



1(8 = 0) - 1(8 = tt) 



(6) 



When the magnetizations of the two FMs are non- 
collinearly arranged, the TMR ratio can be described by 
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FIG. 2: (Color online) Tunnel current as a function of the bias 
voltage (a)-(d) and of temperature (e)-(f) in parallel 1(0) and 
antiparallel I(ir) configurations of magnetizations for differ- 
ent coupling parameter A = 0, 0.1, 0.2, 0.3, 0.5, 0.8, where the 
polarization P = 0.5. 
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Obviously, TMR(n) = TMR, and TMR(0) = 0. 

In the following, for the sake of simplicity for numeri- 
cal calculations, and considering that the electrons near 
the Fermi level in metals are dominant in the tunneling 
process, we may suppose r Q (e) g ' TgT = T ah r a (s) q > lq [ = 
T a i, and the polarization P a = (T a ^ — T a i)/(T a ]- +r a jj. 
If the two ferromagnets are made of the same materials, 

then p L = p r = p, r LT = v m = r, r Li = v Ri = *=f r. 

We will take Iq = eT/h and Go = e 2 /h as scales, respec- 
tively, for the tunnel current and the differential conduc- 
tance, and hereafter take T as an energy scale |2fj|. 

The bias- and temperature-dependence of the tunnel 
current in the parallel and antiparallel configurations of 
magnetizations are presented for different coupling con- 
stant A of the MFL, as shown in Fig. 2. It is seen that 
when A = 0, namely, the MFL recovers to the normal 
Fermi liquid in this case, the tunnel current is propor- 
tional to the bias voltage at small bias, suggesting that 
the system behaves an ohmic law in this case, in agree- 
ment with the conventional result in the Fermi liquid. 
With increasing the coupling constant A, I ~ V curves 
deviate obviously the linear relation, and non-ohmic be- 
haviors appear, i.e. the current increases quadratically 
with the bias voltage. The larger the coupling A, the more 
obvious the distinction from the ohmic behavior, as illus- 
trated in Figs. 2(a)-(d). This observation shows that the 
interactions between electrons in the normal metal would 
have a remarkable effect on the current- voltage character- 
istics where the Ohm law no longer holds. An alternative 
reason for the nonlinearity of / — V characteristics may 




FIG. 3: (Color online) Differential conductance as a function 
of bias voltage in parallel G(0) and antiparallel G(tt) configu- 
rations for different coupling constant A = 0.1, 0.2, 0.3, 0.5, 0.8 
at temperature T = (a)-(b) and T = 0.2Y/k B (c)-(d), where 
the parameters are taken the same as in Fig. 1. 



be that the energy dependent self-energy of the MFL in 
the central region leads to a renormalization of the den- 
sity of states which becomes energy dependent, thereby 
resulting in a nonlinear voltage dependence of the cur- 
rent. When A is small, the tunnel current almost does 
not change with temperature; while A becomes larger, 
the current increases slowly with temperature, as shown 
in Figs. 2(e)-(f). This behavior also differs from that 
in the usual Fermi liquid where the current decreases 
slowly with increasing temperature, as thermal fluctu- 
ations enhance scatterings of conduction electrons and 
thereby contribute to the resistance of the system. It is 
interesting to note that the typical I—V characteristics of 
Ni 80 Fe 2 o/Co/Al-oxide junction (Figure 3.10 in Ref.^l 
are very similar to the shapes of the curves shown in our 
Figs. 2(a)-(d). 

The differential conductance can be obtained by G = 
dI(V)/dV. The results are shown in Figs. 3(a)-(d). As 
A = 0, the conductance is independent of the bias voltage, 
which is nothing but the Ohm law. When A is nonzero, 
the differential conductance behaves as G = Go + G\V 
with Go and G\ nonzero constants at low biases. The 
non-ohmic behavior of G comes from the interactions be- 
tween conduction electrons via the exchange of charge 
and spin fluctuations in the central region. The differ- 
ential conductance is observed to increase slowly with 
increasing temperature at larger A, and almost does not 
change when A is smaller (e.g. A = 0.1). This observation 
is manifested itself in Figs. 2(e)-(f). We notice that the 
linear bias-dependence of the differential conductance in 
various of junctions with Lai.ssSro.isCuO^In and 
even YBCO films have also been observed. It is wor- 
thy of noting that the differential conductance of a con- 
tact between an ordinary metal and a MFL is shown to 
depend linearly on the applied voltage[24|, where due to 
the asymmetry of electrodes, the conductance for positive 
and negative biases is asymmetric. This result is com- 
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FIG. 4: (Color online) Tunnel magnetoresistance as a function 
of the bias voltage as a function of bias voltage at T — 
(a) and T — 0.2r/fc_g (b) and as a function of temperature 
(c) for different coupling constant A = 0.1,0.2,0.3,0.5,0.8 at 
V — 5F/e, where the parameters are taken the same as in 
Fig. 1. 



patible with our observation. The origin of the linearity 
between the conductance and the bias voltage could be 
explained by assuming charging effects j^, the voltage- 
dependent tunneling penetration probabilities 26], DOS 
effects 0,H3, inelastic scattering^, and so on. Our 
present study might offer a different possibility, namely, 
such a linearity between G and V could result from strong 
interactions between conduction electrons via exchanging 
the charge and spin fluctuations. Since the real part of 
the self-energy gives the correction of the single-particle 
energy, describing the elastic scattering of quasiparticles, 
whereas the imaginary part determines the lifetime of the 
quasiparticles, reflecting the inelastic scatterings. There- 
fore, the linearity between G and V could also be dom- 
inated by the inelastic scatterings between conduction 
electrons. 

The TMR ratio as a function of the bias and temper- 
ature for different coupling constant A is shown in Figs. 
4(a)-(c). It is seen that the TMR decreases with increas- 
ing the absolute magnitude of the bias, and is symmetric 
to the zero-bias axis. The larger the coupling constant 
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FIG. 5: (Color online) Tunnel current as a function of the 
coupling constant A at T = (a) and T = 0.2F/k B (b); Tun- 
nel magnetoresistance as a function of the coupling constant 
A at T = (c) and T = 0.2F/k B (d), where V = 5F/e, and 
the other parameters are taken the same as in Fig. 1. 



A, the more rapidly decreasing the TMR, as presented 
in Figs. 4(a) and (b). It suggests that the strong inter- 
actions between conduction electrons tend to suppress 
the TMR ratio, which is a disadvantage for the applica- 
tion of the FM-MFL-FM tunnel junction as a possible 
MR AM. This property of the TMR has also been ob- 
served in various junctions (see Figure 3.7 in Ref . |2l| ) . 
One may observe that the TMR decreases slowly with 
increasing temperature, as shown in Fig. 4(c). 

The current and the TMR ratio as functions of the 
coupling constant A in the MFL for different tempera- 
tures are presented in Figs. 5(a)-(d). It is found that the 
current depends linearly on the coupling constant A in 
the parallel or antiparallel alignment of magnetizations. 
This behavior is also manifested in Figs. 2(a)-(d). It can 
be understood that, with the increase of the coupling 
constant, the single-particle scattering rate which is pro- 
portional to A, increases, leading to that the quantum 
well levels in the MFL could be broadened. Such a level 
broadening could make more electrons tunnel through 
the barrier, thereby resulting in an increase of the cur- 
rent with A, as observed in Figs. 5(a) and (b). In either 
case, T = or T > 0, 1(0) is greater than I(n), imply- 
ing a spin valve effect (see below). The TMR ratio is 
found to decay with increasing the coupling constant A, 
as shown in Figs. 5(c) and (d), suggesting that the inter- 
actions of electrons are detrimental to the TMR effect. 
This may be that the inelastic scatterings of electrons via 
exchanging the charge and spin fluctuations weaken the 
spin-dependent scattering of electrons, leading to that 
the TMR ratio decreases with increasing A. 

The relative angle 9 dependences of the current as well 
as the TMR ratio for different coupling constant are pre- 
sented in Figs. 6(a)-(d). The current as a function of 
9 shows a cosine-like shape, G(9) ~ Go + G\ cos 9 with 
Go, G\ constants, i.e., it decreases with increasing 9 from 
zero to 7r, as illustrated in Figs. 6(a) and (b) for T = 
and 0.2r/fcs, respectively. The TMR ratio as a function 
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FIG. 6: (Color online) Tunnel current as a function of the 
relative orientation angle 6 at T — (a) and T = 0.2r/fcs 
(b); Tunnel magnetoresistance as a function of the relative 
orientation angle 6 at T — (c) and T — 0.2F/fc_B (d), where 
the coupling constant A = 0,0.1,0.2,0.3,0.5,0.8, V = 5T/e, 
and the other parameters are taken the same as in Fig. 1. 

of 9 shows a shape similar to (1 — cos 9). These results 
display nothing but the spin- valve effect. However, as 
discussed above, the coupling constant A tends to sup- 
press the TMR effect. 

In summary, we have discussed the spin-dependent 
transport in FM-MFL-FM tunnel junctions. It is found 
that the current- voltage characteristics in this system de- 
viate obviously from the ohmic behavior, and the tunnel 
current increases slightly with temperature, which are in 
contrast to those of the system with a Fermi liquid where 
the Ohm law is satisfied. The TMR is observed to de- 
cay exponentially with increasing the bias voltage, but 
to decay slowly with increasing temperature. These re- 
sults are qualitatively consistent with the experimental 
observations found in various junctions, suggesting that 
the present study might offer a possible different route 
to understand the unusual experimental results of the 
I — V and G — V characteristics. With increasing the 
coupling constant of the MFL, the current is shown to 
increase linearly, while the TMR is seen to decay slowly. 
It appears that the interactions between electrons in the 
central normal metal via exchanging the charge and spin 
fluctuations tend to suppress the TMR effect. In addi- 
tion, the spin-valve effect is also observed. 
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